Detection of primordial non-Gaussianity (/nl) i n the WMAP 3-year data at above 

99.5% confidence 
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We present evidence for the detection of primordial non-Gaussianity of the local type (/nl), using 
the temperature information of the Cosmic Microwave Background (CMB) from the WMAP 3-year 
data. We employ the bispectrum estimator of non-Gaussianity described in [1] which allows us 
to analyze the entirety of the WMAP data without an arbitrary cut-off in angular scale. Using 
the combined information from WMAP's two main science channels up to imax — 750 and the 
conservative KpO foreground mask we find 27 < f*NL < 147 at 95% C.L., with a central value of 
/nl = 87. This corresponds to a rejection of /nl = at more than 99.5% significance. We find that 
this detection is robust to variations in Imax, frequency and masks, and that no known foreground, 
instrument systematic, or secondary anisotropy explains our signal while passing our suite of tests. 
We explore the impact of several analysis choices on the stated significance and find 2.5 a for the 
most conservative view. We conclude that the WMAP 3-year data disfavors canonical single field 
slow-roll inflation. 
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PACS numbers: 

It is now widely accepted that tests of primordial non- 
Gaussianity, parameterized by the non-linearity param- 
eter /nl , promise to be a unique probe of the early 
Universe [2[ beyond the two-point statistics. Although 
the non-Gaussianity from the simplest inflation models 
is very small, / NL ~ 0.01 - 1 UMM, 

there is a very 

large class of more general models, e.g., models with mul- 
tiple scalar fields, features in inflation potential, non- 
adiabatic fluctuations, non-canonical kinetic terms, de- 
viations from the Bunch-Davies vacuum, among others, 
that predict substantially higher level of primordial non- 
Gaussianity (see [6[ for a review and detailed references). 

Recent calculations of the perturbations arising in the 
ekpyrotic or cyclic cosmological scenarios |7[ have con- 
cluded that these scenarios can predict /nl much larger 
than single field slow-roll inflation [8]. Detailed calcu- 
lations in these models are fraught with difficulties con- 
nected to matching the perturbations through the cosmo- 
logical singularity at the bounce. However, the current 
calculations suggests that primordial non-Gaussianity of 
the /nl type could be a powerful discriminant between 
ekpyrotic models and standard slow-roll inflation. As 
such, the search for primordial non-Gaussianity is com- 
plementary to the search for the inflationary gravita- 
tional wave background. We will argue in this letter that 
the WMAP 3-year data already distinguishes /nl = 100 
from /nl ~ at a statistically significant level. 

Primordial non-Gaussianity can be described in terms 
of the 3-point correlation function of Bardeen's curvature 
perturbations, 3>(fc), in Fourier space: 

($(k 1 )(k 2 )(k 3 )> = (2n) 3 S 3 (k 1 + k 2 + k 3 )F(h,k 2 , k 3 ).(l) 

Depending on the shape of the 3-point function, i.e., 
F(ki,k2,ks), non-Gaussianity can be broadly classified 
into two classes @. First, the local, "squeezed," non- 



Gaussianity where F{k\, fe, fc 3 ) is large for the configura- 
tions in which k\ <C &2, &3- Second, the non-local, "equi- 
lateral," non-Gaussianity where F(ki,k2,ks) is large for 
the configuration when k\ ~ fc 2 ~ fc 3 . 

The local form arises from a non-linear relation be- 
tween inflaton and curvature perturbations [3|, |4( , curva- 
ton models , or the ekpyrotic models Q . The equilat- 
eral form arises from non-canonical kinetic terms such as 
the Dirac-Born-Infeld (DBI) action [ll|, the ghost con- 
densation [12], or any other single-field models in which 
the scalar field acquires a low speed of sound [13[ . While 
we focus on the local form in this letter, it is straightfor- 
ward to repeat our analysis for the equilateral form. 

The local form of non-Gaussianity may be 
parametrized in real space as 

$(r) = $ i (r)+/ NL ($|(r)-($|(r)}) (2) 

where /nl characterizes the amplitude of primordial non- 
Gaussianity. Note that the Newtonian potential has the 
opposite sign of Bardeen's curvature perturbation, <£. 

The first fast bispectrum based /nl estimator using 
temperature anisotropics alone was introduced in [18[ 
(the KSW estimator). The idea of adding a linear term 
to reduce excess variance due to noise inhomogeneity 
was introduced in [19[. Applied to a combination of 
the Q, V and W channels of the WMAP 3-year data 
up to imax ~ 400 this estimator has yielded the tightest 
constraint on /nl so far: —36 < /nl < 100 (2a) pH ]. 
This estimator was generalized to utilize both the tem- 
perature and E-polarization information in [l| , where we 
pointed out that the linear term had been incorrectly 
implemented in Eq. 30 of [19[. The corrected estima- 
tor enables us to analyze the entire WMAP data without 
suffering from a blow-up in the variance at high i. 

Our analysis. We assume a standard Lambda CDM 
cosmology with following cosmological parameters: — 
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FIG. 1: We show the measured value of the non- linear cou- 
pling parameter /nl using WMAP 3-year maps, and the cor- 
responding 95% error bars derived from the Gaussian simula- 
tions. For this analysis the WMAP KpO mask was used. The 
analysis is done for 4 combinations of the frequency channels: 
coadded Q+V+W, coadded V+W, V, and W. 

0.042, n cdm = 0.239, n L = 0.719, h 0.73, r = 0.09, 
and n s = 1. We will discuss the effect of varying these 
fiducial parameters below. 

We used the generalized bispectrum estimator of pri- 
mordial non-Gaussianity of local type described in 
The generalized estimator is given by 

q _|_ qlinear 
JNL , (Oj 

where N is the normalization factor and S vrirfl and 

Sprim r are ^ ne so caue d trilinear and linear term of the 
estimator respectively. The trilinear term captures the 
bispectrum information about /nl while the linear term 
has vanishing expectation and is designed to reduce the 
scatter in the trilinear term induced by the foreground 
mask and WMAP's anisotropic scan strategy. 

Although our estimator [1] can utilize both the tem- 
perature and E-polarization information of the CMB to 
constrain primordial non-Gaussianity, we have used only 
temperature information of the WMAP 3-year data. For 
the analysis we used various combinations of 8 channels of 
WMAP 3-year raw data: Ql, Q2, VI, V2, Wl, W2, W3, 
and W4. For all the simulations we used the WMAP 
3-year maps in HEALPix format with N p i x = 3145728 
pixels. We focused on the V and W bands, which are the 
main WMAP CMB science channels suffering least from 
foreground contamination. We also applied our estimator 
to Q and Q+V+W to assess sensitivity to foregrounds. 

We performed Monte Carlo simulations to assess the 
statistical significance and errors of our /nl estimates. 
For example for the Q+V+W coadded simulated map, 
we first simulated 8 Gaussian maps using the noise and 
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TABLE I: Non- linear coupling parameter /nl using the 
V+W, Q, and Q+V+W WMAP 3-year raw maps, as a func- 
tion of maximum multipole used in the analysis £ m ax and 
mask Kpl2, Kp2, KpO, and Kp0+ (corresponding f s k y is 
stated in the text and the masks are shown in Fig 2) . The last 
row (750*) shows the mean /nl estimated from Gaussian sim- 
ulations including the WMAP foreground model. Foreground 
contamination biases /nl negatively by similar amounts in 
both the data and the model. 



beam properties of the corresponding 8 channels. Then a 
single map was obtained by pixelwise averaging of these 
8 maps. The same procedure was followed to obtain sim- 
ulated coadded maps of the other channel combinations. 
The Sab and Sbb weight maps for the linear estimator 
[lH were obtained using 800 Monte Carlo simulations 
that include the WMAP noise and foreground masks. 

Figure [T] shows the measured value of the non-linear 
coupling parameter /nl for 4 combinations of coadded 
frequency channels (Q+V+W, V+W, V, and W) as a 
function of maximum multipole £ m ax used in the analysis. 
All the analyses in this figure use the KpO mask. The 
figure shows the 95% C.L. error bars derived from Monte 
Carlo simulations. 

For the coadded V+W map there is evidence of pri- 
mordial non-Gaussianity at more than 95% C.L. for all 
£max > 450. For the coadded Q+V+W map there is a de- 
tection of primordial non-Gaussianity at more than 95% 
C.L. for all £ ma x > 500. Residual suboptimality of our 
estimator results in a larger error bar for the Q+V+W 
combination compared to the V+W combination. 

Using the coadded V+W (the least foreground con- 
taminated) channel with £ m ax — 750, we find 

27 < f NL < 147 (at 95% C.L.). (4) 

This rules out the null hypothesis of Gaussian primordial 
perturbations at 99.5% significance. 

Our analysis provides the most information to date on 
the primordial non-Gaussianity of the local type. For the 
sake of comparison with the previous best result in the 
literature (-36 < /nl < 100, for the coadded Q+V+W 
map at the 2<r level for £ m ax 

«4oo HE EES). 

our con- 
straints using the coadded Q+V+W map truncated at 
Zmax = 400 are: -20.84 < f NL < 83.4 (at 95% C.L.). 
We may conclude that the additional information uncov- 
ered by the Yadav et al. estimator [1] at £ > 400 is 



important for our result. As calculated by Creminelli et 
al. [20[ and verified in simulation by [2l|, there is a con- 
tribution to the estimator variance due to non-zero /nl- 
This widens the confidence interval of the estimator by 
3%. It does not however modify the significance of our 
rejection of the Gaussian null hypothesis. 

Interpretation. A detection of non-Gaussianity has 
profound implications on our understanding of the early 
Universe. We will now argue based on an extensive suite 
of null tests and theoretical modeling that our results are 
not due to any known systematic error, foregrounds or 
secondary anisotropy. 

Since our estimator is based on three-point correla- 
tions, any mis-specification of the WMAP noise model 
would not bias our estimator, since Gaussian instrument 
noise has a vanishing three-point function. Similarly, if 
the CMB were Gaussian, asymmetric beams cannot cre- 
ate non-Gaussianity. Beam far-side lobes can produce a 
small level of smooth foreground contamination at high 
galactic latitude [22] at £ < 10. This effect has been 
corrected in the 3-year maps [23]. Since our signal is not 
frequency dependent this is clearly not a dominant ef- 
fect. Even so, we checked for this or any other large scale 
anomaly (such as the axis of evil) by deleting modes with 
£ < 20 from our analysis. We find that our estimate in- 
creases to /nl = 135 ± 96 at (95 % C.L.), leaving the 
statistical significance of our signal at a similar level. 

To test for non-Gaussian time-domain systematics or 
non- Gaussian noise, we take difference between the pairs 
of yearly WMAP data. This creates 3 jackknife realiza- 
tions of WMAP noise maps for each detector including 
real instrument systematics. Applying our estimator to 
these maps gave negligible /nl (~ 1) for all three pairs 
of years, leading us to a conservative bound on the sys- 
tematic error arising from such effects of ±2. 

Seeing the same behavior in all channel combinations 
suggests that even if the detection of non-Gaussianity 
were not primordial, its source would not be frequency 
dependent and not associated with the main galactic 
foregrounds near the galactic plane. This is confirmed 
by repeating our analyses with several foreground masks 
(WMAP Kpl2 mask with f sky = 94.2%, WMAP Kp2 
mask with f sky 84.7%, WMAP KpO mask with f sky 
76.8%, and a larger mask (Kp0+) with f sky = 64.3%, 
which was obtained by smoothing and thresholding the 
KpO mask) for computing the value of /nl using the Q, 
Q+V+W and V+W maps. The results for these analy- 
ses are shown in Table HI It is clear from these analyses 
that the Kpl2 and Kp2 masks do not exclude galactic 
foregrounds at the required level. The Q+V+W combi- 
nation is foreground contaminated even for KpO mask, as 
is also clear from Q band analysis. However, for V+W, 
increasing the mask beyond KpO does not change the re- 
sults significantly. 

To quantify the level at which foreground contamina- 
tion might be expected to affect our results, we per- 
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FIG. 2: The various masks used for computing the /nl in 
Table U. From left to right, and top to bottom we show Kpl2, 
Kp2, KpO, and Kp0+. The point source exclusion regions are 
identical to those in the WMAP. Enlarging the point source 
exclusions does not change our results appreciably. 

form a null test (see last row of Table IJ). We apply 
our analysis to simulated Gaussian CMB skies, expect- 
ing to measure /nl consistent with 0. We model the 
expected diffuse synchrotron, dust and free-free emission 
the same way the WMAP team did to produce their fore- 
ground cleaned maps. We generate simulated WMAP Q, 
V and W maps including this foreground model, Gaus- 
sian CMB and WMAP noise, and we found negligible /nl 
(— 6 ± 4 standard error in the mean for 70 simulations) 
with imax — 750 and KpO mask. The same analysis run 
on simulated Q+V+W data gives an /nl ~ —13 ±5. /nl 
becomes increasingly negative for smaller masks, follow- 
ing the same pattern seen in the real data, and is negli- 
gible for the Kp0+ mask. 

Note that in the cases that we expect to be affected by 
foregrounds (e.g. small masks), /nl is biased negatively 
both in the real data and in simulations. It is therefore 
plausible that any bias due to residual foreground con- 
tamination would falsely reduce the significance of our 
detection. We assign a systematic error of ±[j t° diffuse 
foreground effects in V+W, kpO and £ ma x = 750. 

To get an additional handle on foregrounds, we have 
also analyzed the WMAP foreground-reduced maps. Ir- 
respective of the mask used we always find that fore- 
ground subtraction increases the /nl estimate. How- 
ever, since we cannot guarantee that cleaning the fore- 
grounds does not oversubtract foregrounds we conserva- 
tively quote the results from the raw maps. 

For both the KpO and the extended masks, increasing 
imax increases the central value of /nl- This effect may 
also be seen in Fig. 1. This change is somewhat larger 
than what we measure in Monte Carlo simulations, but 
attaching a significance to this observation is complicated 
by its a posteriori nature. In any case, foregrounds are 
likely not responsible since most models of the high lati- 
tude galactic foreground emission predict more rapid de- 
cay of spatial structure in these foregrounds compared 
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to the CMB on the relevant scales. Note that the same 
arguments also apply to a spinning dust component. 

Our added statistical power relies on information at 
relatively small angular scales, where unresolved point 
source emission may contaminate the map. There are 
several reasons why it is unlikely to be responsible for 
our results. First, jointly fitting for both the primordial 
and the point source bispectrum contributions has shown 
that neglecting the point source contribution does not 
bias /nl on the relevant range of scales 0, [H, Q . Sec- 
ond, since bright point sources are highly clustered, one 
could hypothesize that there may be a contribution due 
to clustered point sources that are not explicitly masked 
by the KpO mask but that are near the point source ex- 
clusion regions. A re-analysis of the data using a mask 
with larger point source exclusion regions showed no sig- 
nificant shift. Galactic point sources with a gradient in 
galactic latitude cannot be responsible because our esti- 
mate is insensitive to extending the mask. 

All studies of the dominant secondary anisotropics con- 
clude that they are negligible for the analysis of the 
WMAP data for l max < 800 @, The largest ex- 

pected biases arise from the ISW-lensing and SZ-lensing 
contributions. These partially cancel due to their oppo- 
site signs and are expected to contribute a net bias of 
~ 2, and a much smaller effect on the error bar. Masking 
the Cold Spot found by [25[ we find the central value of 
/nl increases to (~ 94), but we ignore this enhancement 
of our signal as an a posteriori effect. 

Finally we study the dependence of /nl on cosmologi- 
cal parameters. We find that within the 2a allowed range 
of cosmological parameters, the central value and the er- 
ror bar varies by than 10% with the largest effect due 
to variations in n s . Setting n s = 0.95 reduces all kpO 
estimates by 5 — 20% while simultaneously reducing the 
variance by a similar amount. For n s = 0.95, l max = 750 
has a significantly smaller variance than those for lower 
lmax and gives /nl = 83.5 ± 27, increasing the statistical 
significance of the detection to more than 3a. 

We conclude that the WMAP 3 year data contains 
evidence that allows us to reject the null- hypothesis of 
primordial Gaussianity at the 99.5% significance level. 
Including our systematic error estimates, our result dif- 
fers from = at the 2.5<r level. If our result holds up 
under scrutiny and the statistical weight of future data, 
it will have profound implications on our understanding 
of the physics of the early Universe. As it stands, the 
data disfavors canonical single field slow-roll inflation. 

In addition to repeating our analysis on future data, 
further tests on currently available data using different 
higher order moments and additional bispectrum config- 
urations may provide additional clues. The implications 
of our result for other probes of the cosmological density 
field, such as the mass function and large scale structure 
data should also be considered. This detection demon- 
strates the promise of targeted searches for primordial 
non-Gaussianity as a probe of the early Universe. 
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